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摘 要 : 目前 遥感 监测 土壤 含水 率 的 方法 较 多 ,本 次 研究 选取 了 稳定 性 较 好 、 应 用 较为 广泛 且 所 需 
气象 资料 少 的 3 种 监测 方法 一 一 SWEPDI 指数 法 、 能 量 指数 法 与 TVDI 指数 法 。 以 2016 年 4 月 和 9 
月 两 期 中 高 分 辩 Landsat8 数据 为 数据 源 ,分 别 将 SWEPDI 光谱 法 能量 指数 法 与 TVDI 指数 法 按 不 
同时 间 、 不 同 土 层 深度 与 对 应 时 间 的 土壤 含水 率 野 外 实测 数据 进行 线性 拟 合 ,并 进行 各 模型 之 间 的 
比较 ,选择 出 更 加 适当 的 模型 。 同 时 利用 景观 指数 在 斑 块 类 型 上 分 析 2 ws 
结果 表明 :7TVDI 法 效果 明显 优 于 SWEPDI 光谱 法 和 能 量 指数 法 ,同时 该 方法 还 解决 了 其 他 方法 不 
能 连续 监测 土壤 含水 量 的 问题 ,适用 于 各 种 植被 覆盖 条 件 下 以 及 各 种 土 层 深度 的 土壤 水 分 反 演 。 
此 外 ,选用 6 个 景观 指数 分 析 了 2016 年 4 月 与 9 月 不 同 旱情 等 级 干旱 等 级 的 景观 格局 变化 ,发 现 4 
月 份 PLAND 指数 达到 了 58.76% ,而 9 月 份 轻 早 等 级 的 PLAND 指数 达到 了 44.16% ,都 占据 优势 
地 位 ,其 中 LPI AREA_CV AI 指数 的 值 也 都 达到 了 最 大 ,其 旱情 状况 有 了 明显 的 改善 。 
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土壤 水 分 是 反映 荒漠 化 土地 环境 演变 的 重要 ”针对 低 植被 覆盖 的 区 域 ; 蒸 散 法 主要 利用 气象 资料 
因子 ,是 监测 土地 退化 的 重要 指标 ,也 是 气候 环 ”协同 作物 冠 层 温 度 反 演 土 壤 含 水 率 ,同样 涉及 大 量 
境 . 生 态 状 况 以 及 农 牧 业 生产 等 领域 的 主要 参 参数 ,实现 过 程 复杂 ;植被 指数 法 借助 于 植被 指数 与 
"TOU ,决定 着 区 域 荒漠 化 土地 的 发 生 ` 演 变 以 及 农 ”植被 覆盖 下 的 土壤 含水 率 相关 性 实现 土壤 水 分 监 
牧 业 生产 力 等 方面 2 。 由 于 土壤 含水 率 在 时 空 分 。” 测 , 有 一 定 的 滞后 性 ,不 利于 及 时 获取 田间 信息 。 利 
布 情况 上 变化 极 大 ,传统 的 土壤 含水 量 监测 方法 用 植被 指数 与 温度 结合 “进行 土壤 湿度 的 反 演 ,是 
野外 采集 数据 慢 .消耗 人 力 物 力 多 、 仅 限于 单 点 代 ”对 于 土壤 干 湿 状 态 的 一 种 新 理解 , 现 已 被 广泛 研究 
表 性 差 等 缺点 ,难以 满足 大 范围 实时 监测 的 需求 ， “应 用 5 。 
从 而 限制 了 它 的 应 用 范围 。 遥 感 方法 监测 土壤 含 遥感 .地 理 信 息 系统 与 景观 生态 学 理论 共同 形 
水 量 的 优势 在 于 能 以 不 同时 空 尺 度 不 断 地 观测 地 ”成 了 对 大 尺度 生态 系统 空间 格局 进行 研究 的 独 具 特 
表 状 况 ,提供 地 表 特 征 信 息 , 具 有 监测 范围 广 . 速 色 的 研究 模式 ,利用 景观 指数 对 空间 格局 进行 定量 
度 快 .成 本 低 , 能 够 长 期 对 研究 区 域 的 地 表 特 征 信 ”化 研究 是 这 种 研究 模式 的 基本 内 容 。 国 内 外 的 许多 


息 进行 动态 观测 。 学 者 已 利用 景观 指数 对 不 同 地 区 ` 流 域 的 植被 景观 
利用 遥感 手段 监测 土壤 湿度 的 方法 有 热 惯量 《格局 进行 了 分 析 ,并 取得 了 良好 的 效果 “ 
法 ZIOA .植被 指数 法 .SWEPDI 指数 法 、 能 量 指数 乌 审 旗 地 处 毛乌素 沙 地 腹部 。 近年 来， 干旱 


法 ,温度 植被 指数 法 ” 。 热 惯量 法 通过 求解 地 表 热 。 灾害 发 生 频 率 越 来 越 高 , 牧 革 和 畜产 品 减少 ;草地 
传导 方程 来 实现 ,表达 形式 多 样 . 参 数 复杂 , E ER 退化 和 荒漠 化 加 剧 ; 人 冀 饮 水 困难 ,生命 遭 到 威 


© 收 稿 日 期 : 2018 -04 -12; 修订 日 期 : 2018 -07 -24 
基金 项 目 : 国家 自然 科学 基金 (51769021 ,51169016 ) ;内 蒙古 自然 科学 基金 (2015MS0513) 
作者 简介 : 王 思 楠 (1993 - ) , 男 , 汉 族 , 山 西 朔 州 人 ,硕士 研究 生 ,主要 从 事 水 利和 遥感 信息 技术 研究 . E-mail :nmgnydx2016@ 163. com 
通讯 作者 : 李 瑞 平 (1973 - ) , 男 , 汉 族 ,呼和浩特 人 ,博士 ,教授 ,博士 生 导 师 , 主 要 从 事 节 水 灌溉 与 农业 水 利信 息 技术 研究 . 
E-mail : nmglrp? 163. com 


201810.00134v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


王 思 楠 等 :基于 遥感 数据 对 毛乌素 沙 地 腹部 旱情 等 级 的 景观 变化 特征 分 析 


胁 。 作 为 生态 脆弱 区 ,是 我 国 北方 半 干 旱 区 域 生 
态 变化 研究 的 重点 区 域 ,具有 很 强 的 代表 性 。 
此 ,本 文 以 乌 审 旗 地 处 毛乌素 沙 地 腹部 为 例 ,比较 
分 析 SWEPDI 指数 法 能量 指数 法 与 TVDI 指数 法 
在 乌 审 旗 的 适用 性 ,并 以 实测 的 含水 率 对 所 挑选 
模型 进行 检验 。 以 此 为 基础 ,对 毛乌素 沙 地 ” 景 
观 格局 变化 进行 分 析 , 以 期 为 当地 生态 环境 治理 
提供 科学 参考 。 


1 研究 区 概况 与 研究 数据 


1.1 研究 区 概况 

研究 区 乌 审 旗 处 于 蒙 、. 宁 、 陕 的 三 角 地 带 , 位 于 
内 蒙古 自治 区 最 南端 ,地 势 由 西北 向 东南 倾斜 ,海拔 
一 般 在 1 300 ~ 1 400 m, 108°17' ~ 109? 40' E, 
37°38 ' -39?23'N ,面积 11 645 km? 。 其 所 处 位 置 冬 
冷 夏 热 , 慎 夜 温差 大 ,年 平均 气温 6.8 % 。 受 东南 季 
内 影响 ,降水 时 段 大 多 数 集中 在 7、8.9 H ,年 降水 量 
350 ~400 mm。 多 年 平均 蒸发 量 2 443 mm 的 干旱 
区 ,是 中 国 中 、 东 部 平原 重要 的 生态 屏障 。 沙 万 是 本 
地 区 沙 生 植被 主体 ,该 地 区 集 荒 漠 和 草原 为 一 体 , 土 
地 土质 沙 性 大 趋向 沙漠 化 , 盐 碱 士 使 典型 草原 稀 玻 ， 
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生态 环境 十 分 脆弱 。 
1.2 数据 源 与 处 理 

研究 用 的 影像 数据 由 美国 地 质 调查 局 ( USGS) 
获取 2016 年 4 月 21 日 和 2016 年 9 月 28 日 2 期 履 
盖 乌 审 旗 的 Landsat8 OLL/TIRS 影像 数据 。 两 期 影 
像 都 选择 在 与 野外 数据 同步 采集 的 晴空 无 云 或 者 少 
云 成 像 质量 相对 好 的 天 气 情况 。 同 时 利用 ENVI 软 
件 对 Landsat8 OLI/TIRS 图 像 进 行 辐射 定 标 、 利 用 
Modtran 4 模型 进行 大 气 校正 .几何 精 纠正 等 一 系列 
的 预 处 理 。 
1.3 ”野外 样 点 布设 与 土壤 含水 率 数据 采集 

为 了 使 野外 采样 点 可 以 充分 代表 乌 审 旗 区 域 ， 
结合 研究 区 行政 区 域 空间 分 布 及 道路 交通 状况 ,在 
24 个 样 区 布设 了 120 个 样 点 ,每 个 样 区 (10 km x 10 
km) 中 各 子 样 点 之 间 的 间距 大 致 为 1 km。 考 虑 到 
不 同 深度 的 土壤 含水 率 分 布 状况 ,每 隔 10 cm 深度 
利用 土 钼 对 土壤 含水 率 进 行 一 次 采样 , 共 3 组 重复 
的 0~10 cm 、10 ~20 cm 20 ~30 em 土壤 样品 ,然后 
用 烘 干 称 重 法 测定 其 土壤 含水 量 。 由 于 研究 区 范围 
大 ,野外 试验 采样 于 2016 年 4.9 月 影像 成 像 时 间 的 
前 后 3 d 内 完成 样 点 的 采集 。 
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Fig. 1 Geographical location in the study area and the sampling point distribution 


Ayd 


f r" 
Jí 


" 
J“ 


1 


0 


| 


A 
U 
„U 


nemo 


- 
um 
ral 


图 2 子 样 点 采样 设计 


Fig.2 The sample point sampling design 


2 研究 方法 


遥感 监测 土壤 含水 率 的 方法 较 多 ,本 次 研究 选 
取 了 稳定 性 较 好 、 应 用 较为 广泛 且 所 需 气 象 资料 少 
的 三 种 监测 方法 一 一 SWEPDI 指数 法 .能量 指数 法 
与 TVDI 指数 法 进行 对 比 、 分 析 , 分 别 得 出 3 种 旱情 
指数 的 优 劣 及 其 适用 范围 。 

以 2016 年 4 月 与 9 月 两 期 中 高 分 辨 Landsat8 
数据 为 数据 源 ,分 别 将 SWEPDI 光谱 法 ,能 量 指数 法 
与 TVDI 指数 法 按 不 同时 间 不同 土 层 深度 与 对 应 
时 间 的 土壤 含水 率 野 外 实测 数据 进行 线性 、 指 数 及 
对 数 拟 合 ,并 进行 各 模型 之 间 的 比较 ,选择 出 更 加 适 
当 的 模型 。 以 该 方法 为 基础 ,提取 毛乌素 沙漠 腹部 
3 个 年 份 的 土壤 水 分 并 进行 分 级 ,然后 选取 景观 指 
数 对 分 级 的 影响 进行 景观 格局 计算 和 分 析 。 

2.1 SWEPDI 模型 

(pi Pu) 表示 NIR-Red 特征 空间 上 的 任意 一 
点 ,由 该 点 到 直线 工 的 距离 的 大 小 可 以 说 明 地 表 的 
干旱 情况 , 即 离 L 线 越 远 地 表 越 干旱 ,反之 亦 然 。 
在 NIR-Red 特征 空间 上 任意 取 一 点 已 (pw pw ) 到 
直线 L 的 距离 的 EF 可 以 描述 干旱 的 状况 , 即 距离 
等 于 垂直 干旱 指数 (PD7) 大 小 ,因此 可 以 建立 一 个 
基于 NIR-Red 光谱 空间 特征 的 干旱 监测 模型 PDI 
(图 2) 。 本 文中 采用 郭 兵 修正 后 的 EPDI Bio, 
公式 如 下 : 

Epp] -2 + Mp,, - Lx (DVI - Mp.) (1) 
EVI x s1 +M 
式 中 :pw wu 分 别 为 经 过 FLAASH 大 气 校正 的 红 波 
段 . 近 红外 波段 反射 率 ;WM 为 土壤 线 斜 率 ;EVI 为 增 
强 型 植被 指数 ;L 为 调整 系数 ,按照 研究 区 植被 覆盖 
状况 取 值 ,在 此 取 值 为 0.6。 
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在 土壤 极度 干旱 时 ,EPDI 指数 均 趋向 于 0, 而 
当地 表 湿 度 很 大 或 者 地 表 为 水 体 时 ,EPDI 指数 则 趋 
近 于 1, 因 此 用 1 减 去 各 种 归 一 化 干旱 指数 便 得 到 
了 土壤 水 分 指数 SWEPD , 即 为 : 

SWEPDI =1 - EPDI (2) 
2.2 ”能 量 指数 法 模型 

根据 土壤 热力 学 理论 ,地 球 表面 接收 太阳 辐射 
以 后 ,经 过 复杂 的 转换 又 以 长 波 辐射 的 形式 向 外 发 
射 能 量 ,因此 , 较 强 的 长 波 辐射 标志 着 地 表 温 度 较 
高 。 能 量 指数 的 基本 思路 是 地 球 表面 单位 面积 上 得 
到 的 短波 辐射 为 (1 -A ) ,土壤 越 干燥 ,经 过 转换 向 
外 放出 的 长 波 辐射 越 强 ,表现 为 地 表 和 植 冠 温度 越 
高 ,因为 土壤 中 或 植被 中 的 水 分 吸收 了 一 部 分 太阳 
辐射 ,土壤 越 湿 润 ,经 过 转换 向 外 放出 的 长 波 辐射 越 
弱 ,表现 为 地 表 和 植 冠 温度 越 低 。 利 用 般 感 数据 与 
土壤 增 情 数据 建立 能 量 指数 关系 ,公式 如 下 : 

Dz-z(1-A,)7T; (3) 
式 中 :4 为 近 红外 波段 (光谱 范围 841 ~ 876 nm) 的 
反照 率 ;7, 为 地 表 温 度 。 
2.3 TVDI 模型 

PRICE 46? 人 研究 发 现 植被 指数 和 地 表 温 度 
呈 显 著 的 负 相 关 关 系 ,土壤 水 分 条 件 和 植被 覆盖 情 
况 变化 较 大 时 ,以 遥感 数据 反 演出 的 植被 指数 和 地 
RREN X Y 坐标 的 散 点 图 呈 三 角形 。 在 植被 指 
数 和 地 表 温度 的 三 角形 特征 空间 中 ,可 以 提取 到 干 、 
湿 边 方程 : 


Smax =01 +b; x NDVI (4) 
Smin = Q5 F b, x NDVI (5) 


式 中 :Ts 为 干 边 ,Ts 为 湿 边 ,分 别 由 植被 指数 和 


0.50 


红 光 反射 率 


图 3 NIR-Red 特征 空间 
Fig.3 NIR-Red feature space 


地 表 温 度 根据 干 边 , 湿 边 拟 合 得 到 ;a 、b 、a,、b, 是 
T 湿 边 拟 合 方程 系数 。 

TVDI 依靠 图 像 数 据 由 植被 指数 和 地 表 温 度 
计算 得 到 ,其 定义 为 : 


Smin 


TVDI - 25-1. 
TT 


Smax Smin 


(6) 


TVDI 的 取 值 范围 为 0 ~ 1 ,与 土壤 表层 含水 量 
呈 负 相关 ,与 干旱 程度 呈正 相关 。 即 :TYZD/ 值 越 大 ， 
土壤 湿度 越 低 ,相对 干旱 程度 越 严 重 ; 相 反 , 相 对 干 
FFERR 。 
2.4 景观 格局 分 析 

基于 景观 指标 定量 分 析 景 观 格局 的 特征 与 变化 
是 景观 生态 学 研究 的 核心 之 一 ,合适 的 景观 指数 对 
景观 格局 分 析 的 合理 性 十 分 重 性 。 所 以 本 文 在 借鉴 
前 人 研究 结果 的 基础 上 ,结合 各 指标 生态 学 涵义 和 
研究 区 特点 ,利用 Fragstats 4. 2 软件 毛乌素 沙 地 腹 
部 2 个 时 相 的 景观 指数 ,在 类 型 水 平 上 面积 (C4) 、 
斑 块 数量 (NP) ke E (PLAND) RARER% (LPI) 
斑 块 面积 变异 指数 (4RE4_CY) BERRA BE CAD) 等 
6 个 指标 ,各 景观 格局 指标 的 生态 学 意义 与 计算 公 
式 参见 文献 。 


3 结果 与 分 析 


模型 的 回归 拟 合 
采用 2016 年 4 月 和 9 月 的 2 期 各 土 层 深度 (0 
~10 cm 10 ~20 cm,20 ~30 cm) 的 土壤 含水 率 数 
据 , 分 别 于 SWEPDI 指数 法 能量 指数 法 LTVDI 指数 
法 土 层 深度 含水 率 建立 回归 方程 (WW=a + 80), 回 
归 检 验 结果 见 表 1.2 3。 

总 体 来 看 ,不 同时 期 的 SWEPDI 指数 法 能量 指 
数 法 各 呈现 不 同 程度 的 正 相 关 , 即 指数 越 低 ,土壤 含 
KERR, TVDI 指数 法 呈现 不 同 程度 的 负 相关 , 即 

睛 数 越 高 ,土壤 含水 量 越 低 ,旱情 程度 越 严重 ,因此 
可 以 用 于 表层 土壤 含水 量 和 旱情 监测 研究 。 

4 H fy SWEPDI 指数 法 、 能 量 指数 法 与 土壤 含 
水 率 (0 ~10 cm 10 ~20 em 20 ~30 cm) 的 相关 系数 
R (& ^y i] Jy 0. 606,0. 541,0. 502 ; 0. 652,0. 621, 
0.566 ,整体 上 低 于 9 月 份 的 相关 系数 尺 值 0.641、 
0.594 0. 571,0. 749,0. 670,0. 661, mi 4 月 份 TVDI 
指数 法 与 土壤 含水 率 (0 ~ 10 cm 10 ~20 cm,20 - 30 
cm) 的 相关 系数 R 值 分 别 为 - 0.865、- 0. 840, 
-0.723 整 体 上 高 于 于 9 月 份 的 -0.752 、-0.703、 
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表 1 SWEPDI 对 土壤 含水 率 的 拟 合 结果 
Tab.1 SWEPDI's fitting and test results of soil 


moisture content 


ITA JA. 相关 相对 

二 全 ER B mo 

时 间 ”深度 Bo 系数 误差 
模型 项 o 

/cm n R A% 

0~10 80 线性 0.606 35.381 -17.444 14.96 

4H 10-20 80 ”线性 0.541 32.824 -15.157 21.96 

20 ~30 80 ”线性 0.502 30.017 -12.531 23.45 

0 ~10 80 线性 0.641 50.779 -24.028 14.91 

9H 10-20 80 线性 0.594 41.462 -17.353 18.37 

20 ~30 80 线性 0.571 38.016 -15.094 19.14 


表 2 能量 指数 法 对 土壤 含水 率 的 检验 结果 
Tab.2 The fitting and test results of soil moisture 


content by energy index method 


土屋 样本 二 关 x dex X 

um xx a DU xg 4 "Rx 
/ cm n R i ` A% 

0~10 80 ”线性 0.652 106.470 -18.623 13.84 

4 月 10~20 80 线性 0.621 96.066 -15.523 16.27 
20 ~30 80 线性 0.566 92.419 -14.024 21.41 

0 ~10 80 ”线性 0.749 119.800 -31.108 12.24 

9H 10-20 80 ”线性 0.670 101.620 -24.465 15.67 
20 ~30 80 线性 0.661 92.188 -21.271 16.65 


表 3 TVDI 对 土壤 含水 率 的 拟 合 结果 
Tab.3 The fitting and test results of TVDI on 


soil moisture content 


um xx m PU xa ARR pa 
Zem n 79 R i R A% 

0~10 80 ”线性 -0.865 -26.829 20.202 8.97 

4 月 10~20 80 ”线性 -0.840 -23.312 18.834 10.95 
20 ~30 80 ”线性 -0.723 -21.511 18.399 15.13 

0 ~10 80 ”线性 -0.752 -37.847 27.436 10.85 

9H 10-20 80 ”线性 -0.703 -31.592 24.730 11.01 
20 ~30 80 ”线性 -0.679 -30.848 24.638 14.23 


-0.679, 不管 4 月 份 还 是 9 H TH2S ZRRCR (ED 
随 着 土 层 深度 (0 ~ 10 cm .10 ~ 20 em,20 ~30 em) 的 
增加 而 整体 呈现 出 下 降 趋势 ,说 明 SWEPDI 指数 法 、 
能 量 指数 法 与 TVDI 指数 法 都 适合 0 ~ 10 cm REE 
壤 表 层 含水 率 的 监测 ,这 与 间 峰 等 “认为 Nir-Red 
监测 土壤 水 分 有 效 深度 为 0 ~ 10cm 极为 相符 。 
3.2 验证 

为 了 验证 模型 的 结果 ,在 同期 野外 实测 数据 中 
随机 抽取 40 个 采样 点 的 土壤 含水 量 进行 统计 分 析 
取 和 平均 相对 误差 均 方 根 误差 为 评价 指标 来 进行 定 
量 分 析 验 证 精度 。 相 对 误差 的 平均 值 指标 反映 了 相 


Tos 
m4 验证 结果 
Tab.4 The verification results 
时 间 监测 模型 平均 相对 误差 均 方 根 误差 
SWEPDI 14.31 4.65 
4 月 能 量 指数 法 13.65 4.13 
TVDI 9.68 2.38 
SWEPDI 14.08 4.43 
9 H 能 量 指数 法 13.17 3.98 
TVDI 10.24 2.67 


对 误差 的 整体 状况 , 均 方 根 误差 表示 了 反 演 结果 偏 
离 其 平均 值 的 程度 , 均 方 根 误差 越 大 ,表示 偏离 平均 
值 越 大 。 

从 表 4 验证 结果 可 看 出 :在 4 月 SWEPDI 指数 
法 .能量 指 数 法 .TYD7 指数 法 实测 值 与 反 演 值 的 平 
均 相 对 误差 从 小 到 大 依次 为 9. 68 .13. 65 、14. 31; 均 
方 根 误差 从 小 到 大 依次 为 2. 38,4. 13,4. 65。 在 9 
月 SWEPDI 指数 法 、 能 量 指数 法 、7TVDI 指数 法 实测 
值 与 反 演 值 的 平均 相对 误差 从 小 到 大 依次 为 10. 24、 
13. 17,14. 08; 均 方 根 误差 从 小 到 大 依次 为 2. 67、 
3. 98 4. 43 , 即 平均 相对 误差 . 均 方 根 误差 越 小 稳定 
性 越 高 。 因 此 4 月 .9 月 反 演 结果 稳定 性 从 大 到 小 
依次 为 TVDI 指数 法 .能量 指 数 法 与 SWEPDI 指数 
法 。 
3.3 干旱 指数 适用 性 的 比较 与 分 析 

为 进一步 比较 前 面 分 析 的 三 种 干旱 指数 的 特点 
与 适用 性 ,整理 表 1、2、 和 3 可 以 得 到 表 5。 

由 表 3 5 可 知 : 

(1) 4 H,TVDI 方法 在 0~10 cm 、10 ~20 em, 
20 ~30 cm 深度 所 构建 模型 的 相关 系数 尺 均 大 于 能 


表 5 三 种 方法 的 误差 比较 分 析 
Tab.5 Error analysis of three methods 
e 4 9 BIS 
方法 深度 
SWEPDI 0 - 10 cm 14.96 14.91 14.94 
10 ~20 em 21.96 18.37 20.17 
20 ~30 em 23.45 19.14 21.30 
各 层 平均 A% 20.12 17.47 18.80 
能 量 指 数 法 0~10 cm 13.84 12.24 13.04 
10 ~20 em 16.27 15.67 15.97 
20 ~30 cm 21.41 16.65 19.03 
各 层 平均 A96 17.17 14.85 16.01 
TVDI 0 «10 cm 8.97 10.85 9.91 
10 ~20 em 10.95 11.01 10.98 
20 ~30 cm 15.13 14.23 14.68 
各 层 平均 A96 11.68 12.03 11.86 


注 :表格 中 的 A 为 平均 相对 误差 
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量 指数 和 SWEPDI 方法 。 同 时 ,TVDI 方法 在 0 ~ 10 
cm,10 ~20 cm,20 ~ 30 cm 深度 的 平均 相对 误差 分 
别 为 8.97% 10.95% .15. 1396 ,各 层 平均 相对 误差 
11.68% ;能 量 指数 法 在 0 ~10 cm 10 ~20 cm.、20 ~ 
30 cm 深度 的 平均 相对 误差 分 别 为 13. 84% 、 
16.2796 21. 41% , 各 层 平均 相对 误差 17. 17% ; 
SWEPDI 法 在 0~10 cm 10 ~20 cm .20 ~30 cm 深度 
的 平均 相对 误差 分 别 为 14. 96% 、21. 9696, 
23.45% ,各 层 平均 相对 误差 20. 12% ; 由 此 可 以 知 
道 ,7TVDI 方 法 在 各 层 土壤 含水 率 的 平均 相对 误差 均 
小 于 SWEPDI 法 与 能 量 指数 法 ,并 且 均 在 10% 左右 
拟 合 效果 较 好 。 而 SWEPDI 法 与 能 量 指 数 法 监测 效 
果 则 相对 较 差 。 

(2) 9 H ,TVDI 在 各 层 深 度 所 构建 模型 的 RR 值 
均 比 能 量 指数 法 和 SWEPDI VEKI RRK TV- 
DI 方法 在 0~10 cm.、10 ~20 cm,20 ~30 cm 深度 的 
平均 相对 误差 分 别 为 10. 85% 、18. 01% ,19. 2396 , 
各 层 平均 相对 误差 为 12.03% ;SWEPDI 法 在 0~10 
em,10 ~20 cm,20 ~ 30 cm 深度 的 平均 相对 误差 分 
别 为 14. 91% 、18.37% 、19. 14% ,各 层 平均 相对 误 
差 为 17.47% ;能 量 指数 法 在 0 ~10 cm 10 ~20 em, 
20 ~30 cm 深度 的 平均 相对 误差 分 别 为 12. 2496 、 
15.67% 、16. 65% ,各 层 平均 相对 误差 为 14. 8596 ; 
由 此 可 得 ,7VDI 方 法 在 各 层 的 土壤 含水 率 的 平均 相 
对 误差 均 小 于 SWEPDI 和 能 量 指 数 方法 , 且 各 层 平 
均 相 对 误差 为 12. 03% , 反 演 精度 随 着 土 层 深度 的 
增加 而 降低 ,但 可 以 满足 遥感 反 演 精度 要 求 。 

总 上 所 述 ,3 种 土壤 含水 率 监 测 模型 中 SWEPDI 
模型 的 平均 相对 误差 最 低 而 TVDI 模型 最 高 ,表明 
TVDI 模型 的 土壤 含水 量 估计 值 最 接近 真实 值 而 
SWEPDI 模型 的 偏离 程度 最 大 ;TVDI 模型 的 相关 系 
数 最 高 ,而 SWEPDI 最 低 , 表 明 由 TVDI 模型 反 演 得 
到 的 土壤 含水 量 数据 与 观测 数据 间 的 相关 程度 最 
高 ,而 SWEPDI 最 低 。 同 时 考虑 到 模型 的 简单 .实用 
以 及 稳定 性 ,可 以 确定 TVDI 模 型 是 3 种 模型 中 最 理 
想 的 监测 模型 。 但 是 由 于 能 量 指数 法 计算 公式 中 的 
A, 正 是 植被 对 近 红 外 的 强 反射 区 的 反照 率 ,考虑 到 
了 高 植被 覆盖 时 植被 反射 作用 的 影响 ,同时 也 考虑 
到 地 表 和 植被 冠 层 温度 与 土壤 水 分 成 反比 的 规律 ， 
所 以 能 量 指数 法 也 不 受 应 用 下 垫 面 限制 的 问题 ,对 
土壤 旱情 监测 的 精度 比 SWEPDI 模型 有 所 提高 。 
3.4 ， 斑 块 水 平 景观 格局 分 析 

选取 稳定 性 最 好 的 TVDI 模型 ,提取 毛乌素 沙 
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地 腹部 2016 年 4 月 ,9 月 的 7TVDI 指数 ,并 参照 干旱 
等 级 划分 标准 湿润 (0 < TVDI < 0.2)、 正 常 (0.2 < 
TVDI <0.4) . 轻 旱 (0.4 < TVDI «0. 6) .干旱 (0.6 < 
TVDI «0. 8) 、 重 旱 (0.8 < TVDI <1) ,进行 分 级 (图 
4)" ,从 而 分 析 不 同 旱 情 等 级 的 景观 格局 变化 。 

对 2016 年 4~9 月 斑 块 类 型 上 各 景观 指数 的 变 
化 趋势 进行 分 析 , 由 表 可 看 出 ,C4、NP 从 景观 粒度 
层次 揭示 了 景观 破碎 化 程度 。4 ~9 月 旱情 等 级 出 
现 了 “三 减 两 增 ” 的 现象 ,其 中 湿润 等 级 面积 减少 了 
4 132. 8 ,干旱 等 级 . 重 旱 等 级 面积 减少 了 198 441, 
74 838. 15 ;正常 和 轻 旱 等 级 面积 增加 ,尤其 以 轻 旱 
等 级 的 面积 最 为 显著 ,增加 了 261 405.72, 干旱 等 
级 的 斑 块 数量 NP 从 4 月 的 4 766 个 显著 的 增加 到 
13 705 个 ,增加 了 8 939 个 。 表 明和 干旱 等 级 的 面积 
遭 到 不 同 程度 的 分 割 ,转化 为 低 旱 情 等 级 。 这 一 主 
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要 原因 为 9 月 全 族 降 雨量 稍 多 (图 5 ) 而 导致 ,降雨 
量 为 391 mm, PLAND 度量 斑 块 类 型 占 整个 景观 面 
积 的 相对 比例 ,LPI 是 某 类 型 中 最 大 斑 块 占 整 个 景 
观 面积 的 比例 ,可 以 作为 确定 优势 景观 元 素 的 依据 。 
2016 年 4~9 月 正常 . 轻 旱 等 级 的 PLAND 呈 上 升 趋 
势 , 其 中 轻 旱 等 级 的 PLAND 上 升幅 度 最 大 ,上 升 了 
23.5396 , WHig .干旱 . 重 旱 等 级 的 PLAND $ F 
势 , 其 中 干旱 等 级 的 PLAND 下 降幅 度 最 大 ,下 降 了 
17.87% , 2016 年 4 月 .9 月 各 旱情 等 级 的 LPI 最 大 
值 都 出 现在 PLAND 的 最 大 值 处 。 说 明 低 植被 覆盖 
度 的 4 月 占 优势 旱情 等 级 的 为 干旱 等 级 ;高 植被 覆 
盖 度 的 9 月 占 优势 旱情 等 级 的 为 轻 旱 等 级 。4RE4_ 
CV 表示 同 种 类 型 的 斑 块 面积 的 差异 程度 。4 月 干 
早 等 级 的 斑 块 面积 变异 系数 最 大 ,达到 了 6 009.59。 
9 月 轻 旱 等 级 的 斑 块 面积 变异 系数 最 大 , 达到 了 
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4 2016 年 Landsat8 反 演 TVDI 空间 分 布 状况 
Fig.4 Spatial distribution of TVDI which inversed by Landsat8 in 2016 
表 6 类 型 水 平 下 不 同 旱情 等 级 景观 格局 指数 
Tab.6 The landscape pattern index of different drought conditions under the level of patch 
时 间 等 级 CA PLAND NP LPI AREA, CV AI 
2016.4 湿润 10 140.21 0.9130 1 719 0.031 7 642.417 8 87.060 8 
正常 55 707.12 5.015 1 3 157 0.207 5 860.675 2 90.169 6 
轻 旱 229 137.30 20.628 6 6 179 1.679 4 1 943.3122 93.746 5 
干旱 652 663.62 58.757 7 4 766 25.490 0 6 009. 590 7 97.424 5 
重 旱 163 123.38 14.685 6 4 405 4.7330 4 353.689 5 95.961 3 
2016.9 湿润 6 007.41 0.540 8 2 445 0.073 3 1 575.437 7 91.012 5 
正常 71 713.35 6.456 2 10 792 0.172 9 866. 404 6 87.845 2 
轻 旱 490 543.02 44.162 4 8 964 10. 708 4 5 670.114 6 94.811 1 
干旱 454 222.62 40.892 6 13 705 6.299 3 4 412.444 7 94.672 6 
重 早 88 285.23 7.948 0 11 037 0.618 6 2 184.242 5 91.769 7 
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图 5 2016 年 月 降水 量 分 布 
Fig.5  Rainfall distribution in 2016 


5 670.11。 反 映 出 其 内 部 面积 具有 较 大 的 差异 ,而 
湿润 .正常 .旱情 等 级 的 斑 块 面积 变异 系数 则 一 直 较 
小 ,说 明 其 面积 差异 不 大 。 此 外 ,从 聚合 度 (47) 来 
看 ,4 月 的 干旱 等 级 的 聚合 度 达 97. 42% ,9 月 的 轻 
旱 等 级 的 聚合 度 达 94. 81% ,其 主要 原因 为 这 一 年 
低 植被 覆盖 度 的 4 月 的 降雨 量 为 0 mm, 处 于 干旱 时 
期 (图 5)。 


4 结论 


本 研究 在 SWEPDI 指数 法 .能量 指数 法 与 TVDI 
指数 法 3 种 监测 方法 的 最 优 单一 旱情 指数 选择 的 基 
础 上 ,从 而 分 析 不 同 旱情 等 级 的 景观 格局 变化 。 虽 
然 利 用 7s-NDVI 特征 空间 进行 旱情 监测 已 取得 一 定 
的 进展 ,但 对 于 覆盖 较 好 的 天 然 植被 和 农作物 ,使 用 
NDVI 造成 的 红 光 饱和 问题 比较 严重 ,背景 的 土壤 号 
声 也 在 一 定 程度 上 损害 了 NDVI 的 空间 一 致 性 ,而 
土壤 调节 植被 指数 MSAVI 相对 NDVI 更 能 描述 植被 
覆盖 度 和 土壤 背景 的 优势 。 同 时 考虑 到 毛乌素 沙 地 
腹部 地 形 沙漠 化 .人 为 灌溉 、 降 雨 等 复杂 的 地 理 环 
境 ,为 了 准确 地 反映 干旱 特征 的 目的 ,综合 多 源 数据 
的 干旱 综合 监测 模型 是 研究 复杂 干旱 监测 问题 的 新 
途径 。 本 研究 采用 的 2 个 时 相 的 遥感 资料 , 仅 代表 
着 两 个 不 同时 点 的 现状 水 平 ,如 果 能 够 获得 更 加 连 
续 的 遥感 资料 , 则 对 研究 更 具有 说 服 力 。 

土壤 水 分 胁迫 是 制约 毛乌素 沙 地 生态 环境 和 可 
持续 发 展 的 主要 问题 。 开 展区 域 尺度 下 大 而 积 、 高 
精度 的 土壤 水 分 监测 有 利于 旱情 预报 气象 水 文 等 
领域 研究 。 本 研究 以 毛乌素 沙 地 腹部 乌 审 旗 为 例 ， 
研究 土壤 水 分 的 提取 方法 ,并 以 此 方法 为 基础 ,研究 
沙 地 腹部 土壤 水 分 的 景观 格局 演变 规律 。 得 到 如 下 


(1) SWEPDI 法 与 能 量 指数 法 进行 土壤 含水 率 


监测 的 效果 随 着 时 间 的 变化 而 起 伏 较 大 ,7TVDI 法 监 
测 土壤 含水 率 效 果 的 稳定 性 相对 较 好 。 同 时 ,TVDI 
法 监测 土壤 含水 率 精 度 要 优 于 SWEPDI 法 与 能 量 指 
数 法 ,而 能 量 指数 法 又 稍微 优 于 SWEPDI 法 。 由 此 
来 确定 TVDI 的 线性 模型 为 本 研究 中 最 理想 的 拟 合 

(2) 拟 合 程度 随 着 土 层 深度 (0 ~ 10 cm 10 ~20 
cm .20 ~30 em) 的 增加 而 呈现 出 整体 降低 趋势 , 而 
平均 相对 误差 随 着 土 层 深度 的 增加 而 呈现 出 整体 组 
慢 上 升 趋势 ,说 明 0 ~ 10 em 土 层 深度 的 土壤 含水 率 
监测 精度 比 10 ~20 cm 20 ~30 em 土 层 深度 的 土壤 
含水 率 监测 精度 高 ,但 相差 不 太 大 。 

(3) 4~9 月 旱情 等 级 出 现 了 "三 减 两 增 ” 的 现 
象 ,其 中 湿润 等 级 面积 减少 了 4 132. 8 ,干旱 等 级 . 重 
晶 等 级 面积 减少 了 198 441, 74 838. 15 ;正常 和 轻 旱 
等 级 面积 增加 ,尤其 以 轻 旱 等 级 的 面积 最 为 显著 , 增 
加 了 261 405.72。 其 中 4 月 和 9 月 分 别 以 干旱 等 级 
和 轻 旱 等 级 为 主 , 对 应 的 其 中 ZP1 AREA. CV AI TR 
数 的 值 也 都 达到 了 最 大 。 表 明 相 较 植 被 生长 初期 的 
4 月 而 言 , 在 较 高 植被 覆盖 且 生 长 状况 良好 条 件 下 
的 9 月 具有 更 高 的 土壤 水 分 能 力 。 
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Characteristics of the landscape changes under the hinterland 
drought conditions in Mu Us Sandland 


WANG Si-nan, LI Rui-ping, HAN Gang, WANG Yao-qiang, HU Yong-ping, TIAN Xin 
(Inner Mongolia Agricultural University ,Inner Mongolia , Hohhot 010018 , Neimengu , China ) 


Abstract: At present the methods using remote sensing technology to monitor soil moisture content become popular. 
This paper selected three monitoring methods which are widely applied with better stability and require less meteoro- 
logical data , namely the SWEPDI index method ,the energy index method and the TVDI index method ,to better under- 
stand the characteristics of the landscape changes under the hinterland drought conditions in Mu Us Sandland , Inner 
Mongolia , China. Using Landsat8 data in April 2016 and September 2016 as data source , each of the three methods was 
conducted a linear fitting with the field survey data according to the time ,the soil depth and corresponding soil mois- 
ture content. À comparison of the results from the models determined the more suitable model. At the same time ,the 
change trend of soil moisture in 2 time phases was analyzed based on the plaque type by using the landscape index. 
The results showed that the effect of TVDI index method is superior to both the SWEPDI index method and the energy 
index method ,and the TVDI index method also solved the problem that existed in continuously monitoring the soil 
moisture content. The TVDI index method is suitable for all kinds of vegetation conditions and various soil depth in 
drought monitoring. In addition ,six landscape indexes were chosen to analyze the landscape pattern changes in April 
2016 and September 2016 at different levels of drought ,and it is found that the PLAND reached 58. 7696 in April , and 
it became 44. 16% in September when it is at the light drought grade. Both were dominant as other indexes like LPI, 
AREA. CV and AI index all reached the maximum value. The drought condition was improved obviously. 
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pattern ; 


